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We introduce and motivate the method of effective charges, and consider how to implement an all-orders resummation 
of large kinematical logarithms in this formalism. Fits for QCD A and power corrections are performed for the e+e" 
event shape obesrvables 1-thrust and heavy-jet mass, and somewhat smaller power corrections found than in the usual 
approach employing the "physical scale" choice. 



1. Introduction 

In this talk I will describe some recent work together with Michael Dinsdale concerning the relative size of 
non-perturbative power corrections for QCD event shape observables 0. For e+e" event shape means the DELPHI 
collaboration have found in a recent analysis that, if the next-to-leading order (NLO) perturbative corrections are 
evaluated using the method of effective charges y, then one can obtain excellent fits to data without including 
any power corrections ..S] . In contrast fits based on the use of standard fixed-order perturbation theory in the MS 
scheme with a physical choice of renormalization scale equal to the cm. energy, require additional power corrections 
Ci/Q with Ci ~ 1 GeV. Power corrections of this size are also predicted in a model based on an infrared finite 
coupling Q , which is able to fit the data reasonably well in terms of a single parameter. Given the DELPHI result it 
is interesting to consider how to extend the method of effective charges to event shape distributions rather than means. 



2. The method of effective charges 

Consider an e+e~ observable TZ{Q), e.g. an event shape observable- thrust or heavy-jet mass, Q being the cm. 
energy. 

7^(Q) = aO^ RS) + ri(^/Q, RS)a2(^, RS) -f -frz WQ, RS)a^^i, RS) + • • • , (1) 

Here a = Ofg/Tr. Normalised with the leading coefficient unity, such an observable is called an effective charge. The 
couplant a(/i, RS) satisfies the beta-function equation 

dln[fi) 

Here b = (33 — 2iV/)/6 and c = (153 — 19iV/)/12& are universal, the higher coefficients Ci, i > 2, are RS-dependent 
and may be used to label the scheme, together with dimensional transmutation parameter A [J. The effective charge 
TZ satisfies the equation 

= p{TZ{Q)) = -67e2(l + c7^ + p2n' + p^TZ' + ■■■). (3) 

This corresponds to the beta-function equation in an RS where the higher-order corrections vanish and TZ = a, the 
beta-function coefficients in this scheme are the RS-invariant combinations 

P2 = C2 + r2 - ric - rj 

P3 = C3 + 2r3 - 4rir2 - 2rip2 ~ r\c + 2r\. . (4) 
Eq.(3) for dTZ/d\nQ can be integrated to give 
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Figure 1: Values of A-rj-e obtained for hadronization corrected 1-thrust data. 



The dimensionful constant Atj arises as a constant of integration. It is related to the dimensional transmutation 
parameter by the exact relation, 

At^ = e'^'^'Kjjs ■ (6) 
Here r = ri(l, MS) with /i ~ Q, is the NLO perturbative coefficient. Eq.(5) can be recast in the form 

Aj^=QT{n{Q))g{n{Q))e-^^\2c/bf' . 
The final factor converts to the standard convention for A. Here J' {TV) is the universal function 

and Q{Tl) is 

P2 



(7) 
(8) 



g{n) = 1 - ^7^ + o(7^') 





(9) 



Here p2 is the NNLO ECH RS-invariant. If only a NLO calculation is available, as is the case for e'^e jet observables, 
then Q{TZ) = 1, and 

= QHnQ))e-'''\2c/hf"' . (10) 

Eq.(lO) can be used to convert the measured data for the observable TZ into a value of Aj-fg bin-by-bin. Such an 
analysis was carried out in Ref. 6] for a number of e+e^ event shape observables, including thrust and heavy jet 
mass which we shall focus on here. It was found that the fitted A values exhibited a clear plateau region, away from 
the two-jet region, and the region approaching T = 2/3 where the NLO thrust distribution vanishes. The result for 
1-thrust corrected for hadronization effects is shown in Fig. 1. 

Another way of motivating the effective charge approach is the idea of "complete renormalization group improve- 
ment" (CORGI) 0- One can write the NLO coefhcient ri(/i) as 



ri(^) = 61n-j- 
A' 



Wn- 



MS 
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(11) 
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Hence one can identify scale-dependent /x-logs and RS-invariant "physical" UV Q-logs. Higher coefficients are poly- 
nomials in Ti. 

r2 ^ rl + ric+ (p2 - C2) 

r, = rl + \cr\ + (3p2 - 2c2)ri + (y " y ) ■ (12) 

Given a NLO calculation of ri, parts of r2,r3, . . . are "RG-predictable" . One usually chooses /i = xQ then r\ is Q- 
independent, and so are all the r„. The Q-dependence of TZ{Q) then comes entirely from the RS-dependent coupling 
a{Q). However, if we insist that /i is held constant independent of Q the only Q-dependence resides in the "physical" 
UV Q-logs in ri. Asymptotic freedom then arises only if we resum these Q-logs to all- orders. Given only a NLO 
calculation, and assuming for simplicity that that we have a trivial one loop beta- function /3(a) = —ba^ so that 
a(/i) = l/51n(/z/Ap-g) the RG-predictable terms will be 



n = a{^l) [l + Y,{a{^Ji)r,{^i)A 



(13) 



Summing the geometric progression one obtains 



l_(51n^-61n-^)a(^) 



n{Q) = a(/i)/ 

= l/b\n{Q/kn) , (14) 
The /u-logs "eat themselves" and one arrives at the NLO ECH result TZ{Q) = 1 / b\n{Q / An) . 



As we noted earlier, and as will be discussed by Klaus Hamacher in his talk use of NLO effective charge 
perturbation theory (Renormalization Group invariant (RGI) perturbation theory) leads to excellent fits for e+e~ 
event shape means consistent with zero power corrections, as illustrated in Figure 2. taken from Ref.Q. Given this 
result it would seem worthwhile to extend the effective charge approach to event shape distributions. It is commonly 
stated that the method of effective charges is inapplicable to exclusive quantities which depend on multiple scales. 
However given an observable TZ{Qi, Q2, Q3, ■ ■ • , Qn) depending on n scales it can always be written as 

n = 7^(Ql, Q2/Q1, . . . , Qn/Qi)=n,,x,...xAQi) ■ (i5) 

Here the Xi=Qi/Qi are dimensionless quantities that can be held fixed, allowing the Qi evolution of TZ to be obtained 
as before. In the 2-jet region for e+e~ observables large logarithms L = ln(l/xi) arise and need to be resummed to 
all-orders. 



3. Resumming large logarithms for event shape distributions 

Event shape distributions for thrust (T) or heavy-jet mass {ph) contain large kinematical logarithms, L — \n(l/y), 
where y = (1 - T), p/,., • • •. 

= ALdaL^) + L-'ANLL{aL^) + ■■■ ■ (16) 

a dy 

Here LL, NLL, denote leading logarithms, next-to-leading logarithms, etc. For thrust and heavy-jet mass the 
distributions exponentiate |^ 

fy' 1 da 
Jq <y dy 

+ g2{aTTL) + agsiairL) + + •••) + D{a7r, y) . (17) 
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Figure 2; Fits for Qs(Mz) for means of e^e event shape observables taken from Ref.Q. The quality of the "pure RGI" fits 
on the right is noteworthy. 

Here gi contains the LL and g2 the NLL. C = 1 + 0{a) is independent of y, and D contains terms that vanish as 
2/ — * 0. It is natural to define an effective charge TZ{y') so that 



Ry{y') ^ exp{ro{y')n{y')) . 
This effective charge will have the expansion 

rQ{L)n{L) = ro{L){a + ri{L)a^ + r2{L)a^ + • • •) 
Here ro(L) ~ L^, and the higher coefficients r„(L) have the structure 

— + Lj + • ■ • 



(18) 



(19) 



(20) 



Usually one resums these logarithms to all-orders using the known closed- form expressions for giiaL) and g2{aL), 
where a is taken to be the MS coupling with a "physical" scale choice = Q (MSPS). Instead we want to resum 
logarithms to all-orders in the p{TZ) function (ECH). The form of the pn RS-invariants (Eq.(4)) means that the pn 
have the structure 



Pn — Pn ^ + Pn ^ 



(21) 
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Figure 3: Comparison of the 1-thrust distribution using various NLO approximations in the 2-jet region. The solid curve arises 
from exponentiating the NLO ECH. The dashed curve is obtained by expanding this to NLO in MSPS. The dotted curve is 
an unexponentiated NLO ECH fit. DELPHI data at Q = Mz are plotted. Aj^ = 212 MeV is assumed. 



One can then define all-orders RS-invariant LL and NLL approximations to p{TZ), 

oo 

PllW = -67^2(l + c7^ + ^p^LL"7^") 

n=2 
oo 

PnMH) = -67^2(l + cTl+ J^iPn"^^" + Pn^^L^-^W) ■ (22) 

n=2 

The resummed pnll(72.) can then be used to solve for 7?.nll by inserting it in Eq.(5). Notice that since A-r. involves 
the exact value of ri(l,MS') there is no matching problem as in the standard MSPS approach. The resummed 
Pll{T^) can be straightforwardly numerically computed using 

PLL(x)=/?(a)^ = -6a^, (23) 

with a chosen so that 7^LL(a) = x. The same relation with /3(a) = —ba^{l + co.) suffices for pnll{T^), although in 
this case one needs to remove NNLL terms, e.g. an term which would otherwise be included in p2- This can be 
accomplished by numerically taking limits L ^ oo with LIZ fixed. 

As we have noted a crucial feature of the effective charge approach is that it rcsums to all-orders RG- Predictable 
pieces of the higher-order coefficients, thus the NLO ECH result (assuming c = for simplicity) corresponds to an 
RS-invariant resummation (c.f. Eq.(13).) 

a + rxo" r\a^ + • • • + r^a"+i + • • • . (24) 

Thus even at fixed-order without any resummation of large logs in p{Ti) a partial resummation of large logs is 
automatically performed. Furthermore one might expect that the LL ECH result contains already NLL pieces of the 
standard MSPS result. 

In Figure 3 we show various NLO approximations. Notice that the solid curve, which corresponds to the 
exponentiated NLO ECH result, is a surprisingly good fit even in the 2-jet region, whereas the dashed curve 
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Figure 4: For 1-thrust the ratio of tiie NLL MSPS coefficient at 0(a") "predicted" from the LL ECH result to the exact result 
(diamonds). The triangles show the "prediction" from the NLO ECH result. 

which is the NLO AfS'PS result, has a badly misplaced peak. The all-orders partial resummation of large logs in 
Eq.(15) gives a reasonable 2-jet peak. Figure 4 shows that the NLL MSVS coefficients "predicted" from the LL 
ECH result by re-expanding it in the M5PS couphng are in good agreement with the exact coeffiecients out to 0(a^°). 



4. Fits for KjYg and power corrections 

We now turn to fits simultaneously extracting A-jg-^ and the size of power corrections Ci/Q from the data. To 
facilitate this we use the result that inclusion of power corrections effectively shifts the event shape distributions, 
which can be motivated by considering simple models of hadronization, or through a renormalon analysis |ic| . Thus 
we define 

Rpciv) = RpTiy - Ci/Q) . (25) 

This shifted result is then fitted to the data for 1-thrust and heavy jet mass. e~^e~ data spanning the cm. energy 
range from 44 — 189 GeV was used (see |lj for the complete list of references). The resulting fits for 1-thrust and 
heavy-jet mass are shown in Figures 5. and 6.. 

The ECH fits for thrust and heavy jet mass show great stability going from NLO to LL to NLL, presumably because 
at each stage a partial resummation of higher logs is automatically performed. The power corrections required with 
ECH are somewhat smaller than those found with AfS'PS, but we do not find as dramatic a reduction as DELPHI 
find for the means. This may be because their analysis corrects the data for bottom quark mass effects which we 
have ignored. The fitted value of A-p-g- for ECH is much smaller than that found with AfS'PS, {as{Mz) = 0.106 
(thrust) and 0.109 (heavy-jet mass)). Similarly small values are found with the Dressed Gluon Exponentiation (DGE) 
approach 11]. A problem with the effective charge resummations is that the p(TV) function contains a branch cut 
which limits how far into the 2-jet region one can go. We are limited to 1 — T > O.Q5Mz/Q m the fits we have 
performed. This branch cut mirrors a corresponding branch cut in the resummed gi{aL) function. Similarly as 1 — T 
approaches 1/3 the leading coefficient ro(L) vanishes and the Effective Charge formalism breaks down. We need to 
restrict the fits to 1 — T < 0.18. From the "RG-predictability" arguments we might expect that these difficulties 
would also become apparent for a NNLL Af SPS resummation. One will be able to check this expectation when a 
result for g^^anL) becomes available. 
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Figure 5: Fits to 1-thrust for Ajjg and Ci. Solid 2a error ellipses are for ECH, dashed are MSPS. The arrows show the effect 
of varying the scale between Q/2 < /x < 2Q. 




Figure 6: Fits for heavy-jet mass. 
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Figure 7: The dashed hne corresponds to PM Si , and the soUd hne to the physical scale choice M — jj, 
is PMS2 and is in much better agreement with the data points, [l^ 
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5. Extension to event shape means at HERA 



Event shape means have also been studied in DIS at HERA j^^. For such processes one has a convolution of 
proton pdf's and hard scattering cross-sections, 



Mep ^X,Q)_^r ^^^^^^^ ^^^^Mea ^ X, Q, M) 



dX 



dX 



(26) 



There is no way to directly relate such quantities to effective charges. The DIS cross-sections will depend on a 
factorization scale M , and a rcnormalization scale at NLO. In principle one could identify unphysical scheme- 
dependent \ii{M/ Kj^) and ln(/x/A-jgg), and physical UV Q-logs, and then by all-orders resummation get the M and 
/j,-dependence to "eat itself" . The pattern of logs is far more complicated than the geometrical progression in the 
effective charge case, and a CORGI result for DIS has not been derived so far. Instead one can use the Principle 
of Minimal Sensitivity (PMS) and for an event shape mean {y) look for a stationary saddle point in the (/i, M) 
plane ^sj. It turns that there are large cancellations between the NLO corrections for quark and gluon initiated 
subprocesses. One can distinguish between two approaches, PMSi where one seeks a saddle point in the (/i, M) 
plane for the sum of parton subprocesses, and PMS2 where one introduces two separate scales and /ig and finds 
a saddle point in {^q, fig, M). PMSi gives power corrections fits comparable to MSPS with AI = /i = Q- PMS2 
in contrast gives substantially reduced power corrections. This is shown in Figure 7 for a selection of HERA event 
shape means. Given large cancellations of NLO corrections RG-improvement should be performed separately for the 
q and g-initiated subprocesses, and so PMS2 which indeed fits the data best, is to be preferred. 
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6. Conclusions 

Event shape means in e+e^ annihilation are well-fitted by NLO perturbation theory in the effective charge ap- 
proach, without any power corrections being required. With the usual MSVS approach power corrections Ci/Q 
are required with Ci ~ 1 GeV. Similarly sized power corrections are predicted in the model of Ref. 4]. It would 
be interesting to modify this model so that its perturbative component matched the effective charge prediction, but 
this has not been done. We showed how resummation of large logarithms in the effective charge beta-function p{TV) 
could be carried out for e^e~ event shape distrtibutions. If the distributions are represented by an exponentiated 
effective charge then even at NLO a partial resummation of large logarithms is performed. As shown in Figure 3 
this results in good fits to the 1-thrust distribution, with the peak in the 2-jet region in rough agreement with the 
data. In contrast the MSVS prediction has a badly misplaced peak in the 2-jet region, and is well below the data 
for the realistic value of ^jjg — 212 MeV assumed. We further showed in Figure 4 that the LL ECH result contains 
already a large part of the NLL AfS'PS result. We found unfortunately that p{TV) contains a branch point mirroring 
that in the resummed gi{aL) function. This limited the fit range we could consider. We fitted for power corrections 
and A-pg to the 1-thrust distribution and heavy-jet mass distributions, finding somewhat reduced power corrections 
for the ECH fits compared to A/5PS, with good stability going from NLO to LL to NLL. The suggestion of the 
"RG-predictability" manifested in Figure 4 would be that the NLL ECH result contains a large part of the NNLL 
AfS'PS result. This suggests that the branch point problem which limits the ability to describe the 2-jet peak, would 
also show up given a NNLL analysis. This can be checked once the 53 (aL) function becomes available. Recent work 
on event shape means in DIS was briefly mentioned and seemed to indicate that greatly reduced power corrections 
are found when a correctly optimised PMS approach is used. 
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